Phase-modulated quadrature squeezing in two coupled cavities containing a two-level system
Li Hao-Zhen1, 2, †, Zeng Ran1, Zhou Xue-Fang1, Bi Mei-Hua1, Xu Jing-Ping2, ‡, Yang Ya-Ping2
School of Communication Engineering, Hangzhou Dianzi University, Hangzhou 310018, China
Key Laboratory of Advanced Micro-Structured Materials of Ministry of Education, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China

 

† Corresponding author. E-mail: haozhenli@hdu.edu.cn xx_jj_pp@hotmail.com

Project supported by the Fundamental Research Funds for the Provincial Universities of Zhejiang Province, China (Grant No. GK199900299012-015), the National Natural Science Foundation of China (Grant Nos. 11574229, 11874287, and 11574068), and the Science Fund from the Shanghai Science and Technology Committee, China (Grant No. 18JC1410900).

Abstract

The phase-modulated quadrature squeezing in the system that composed of two coupled cavities interacting with a two-level atom is investigated. The variances of the amplitude and phase quadrature of the output field are calculated. It turns out that the squeezing behaviors of the output field can be obviously modified due to the phase difference of the coupling strengths between the atom and the two cavities. The squeezing in one quadrature (i.e., phase quadrature) can be transferred into another (i.e., amplitude quadrature), or the quadrature squeezing located at the low-frequency region can be transferred into the high-frequency region by modulating the relative phase of the coupling strengths. Furthermore, the effects of the decay mismatch between the two cavities and the coupling mismatch between the atom and the cavities on the quadrature squeezing have been discussed. The results show that both the decay mismatch and the coupling mismatch play a positive role in generating better quadrature squeezing.

PACS: ;03.67.-a;;03.65.Yz;
1. Introduction

Nonclassical light source is an essential element for implementation quantum information processing,[1] which is extensively applied to quantum simulations,[2] quantum networks,[3] quantum communications,[4] and linear optical quantum computations.[5] In recent years, the nonclassical properties of light, such as antibunching,[6] entanglement,[7] and squeezing,[8] have been widely studied in quantum optical systems. Squeezing,[9] which possesses variance of a quadrature amplitude or phase below that of a vacuum or coherent state, is regarded as one of the most peculiar nonclassical phenomena. Such a squeezing plays a crucial role in noiseless communications,[10] gravitational wave detectors,[11] and other highly precision measurements,[12,13] such as precision interferometry and atomic spectroscopy.

The generation of squeezed light was first realized in an atomic vapor based on the four-wave mixing process.[14] After then, a variety of schemes for generating squeezed state have been proposed,[1528] including the optical parametric process,[15,16] optomechanics systems,[1719] time-delayed coherent feedback,[20] periodically modulation,[21] and cavity quantum electrodynamics (QED).[6,2224] The physics behind them is mainly related to the nonlinear optical process. For cavity QED system, a single atom in a suitable setup can work as an effective nonlinear-optical medium operating in the quantum regime. Quadrature squeezing can be observed by trapping a single atom into a weakly driven high-Q cavity. It is shown that the squeezing comes from a cavity-enhanced coherent interaction and an optimal squeezing can be obtained by adjusting the coupling strength.[29,30] However, in these proposals, only the amplitude of the coupling strength is considered, while the influence of the phase effect of the coupling strength on the squeezing properties has been ignored.

It is well known that the phase effect plays a crucial role in quantum mechanics. The physical observation of a quantum superposition state in measurement statistics is much different due to the phase difference between the amplitudes. The physical mechanisms of the phase-controlled quantum feature in an atom–photon coupling system have been investigated.[3136] For example, the single-photon transmission and reflection spectra can be significantly influenced by modulating the phase differences of the coupling strengths between the atom and photons,[35] or between the external driving fields.[36] A single-photon router or switch can be realized based on such phase modulation, which can be used as a building block for various quantum information processing.

Recently, phase-dependent quantum statistical properties in a driven cavity QED system have also been reported.[3739] They indicated that the phase difference of the coupling strengths between the atom and the cavity modes or the relative phase between the external driving fields plays an important role in the generation of quantum light. For instance, the photonic statistical properties of the light field from a cavity that contains two atoms are closely related to the atom’s relative positions (which leads to the phase difference of the coupling strengths between the cavity mode and different atoms). The light field can be tuned from antibunched to bunched or superbunched by modifying the phase difference between the coupling strengths.[37] Similar results have also been shown by coupling two cavity modes to a two-level system[38] or a microresonator optomechanical system.[39] They found that photon antibunching can be optimally modulated by properly engineering the relative phase between the coupling strengths or between the driving fields, which is useful for generating single photon sources or nonclassical fields.

Motivated by these advances, in this paper, the phase-modulated quadrature squeezing of the output field in a system contains two coupled cavities interacting with a two-level atom is investigated. The results show that the squeezing properties depend strongly on the phase and an optimal quadrature squeezing can be obtained by carefully engineering the relative phase of the coupling strengths between the atom and the two cavity modes. The squeezing behaviors in the amplitude and phase quadrature, or located at the low-frequency region and high-frequency region can be transferred into each other by modulating the relative phase between the coupling strengths. Furthermore, this phase-modulated quadrature squeezing also related to the decay mismatch and the coupling mismatch. These mismatch can be applied to generate better quadrature squeezing.

This paper is organized as follows. In Section 2, we introduce the model and the basic theory to describe the squeezing properties of the system. In Section 3, we analyze how the phase difference of the coupling strengths between the atom and the two cavities affects the squeezing properties of the output field. Besides, the influences of the decay mismatch between the two cavities and the coupling mismatch between the atom and the two cavity modes on the phase-dependent quadrature squeezing are explored in detail. Finally, a conclusion is drawn in Section 4.

2. Model and basic theory

The model under consideration is constituted by a two-level atom and two coupled single mode cavities. An external coherent probe field is applied to drive the mode of cavity 1, as shown in Fig. 1. Under the rotation wave approximation, the Hamiltonian of the driven system can be written as

Here, ωa is the transition frequency of the two-level atom, ωc,j is the resonance frequency of the j-th cavity, aj and are the annihilation and creation operators of the j-th cavity mode, respectively. σ = |g⟩ ⟨e| and σz = |e⟩ ⟨e| − |g⟩ ⟨g| are the atomic lowering operator and the z-component Pauli matrix, respectively. gj is the coupling strength between the j-th cavity and the atom. In order to show the influence of phase difference between the coupling strengths on the squeezing properties, gj is assumed be a complex number. J is the coupling strength between cavity 1 and cavity 2, which can be tuned by adjusting the distance of the two cavities. ε and ω are the amplitude and frequency of the input field Sin, respectively.

Fig. 1. Schematic diagram of a two-level atom coupled to two single mode cavities with coupling strengths g1 and g2, respectively. J denotes the coupling strength between cavity 1 and cavity 2, ωa denotes the transition frequency of the two-level atom, Sin and Sout represent the input and output of an external probe field, respectively.

In a frame rotating at the input frequency ω, the Hamiltonian in Eq. (1) takes the form

where Δa = ωaω (Δc,j = ωc,jω) is the detuning between the atom (cavity j) and the input field. The dynamics of the system in the presence of dissipation can be described by Heisenberg equation. The equations of relevant operators can be written as

Here, γa is the atomic decay rate, and κj is the leakage rate of the photons from the j-th cavity to the outside. The cavity operators aj can be obtained directly from Eqs. (3) and (4) in the steady state (),

Using the standard input-output theory,[40,41] the output Sout from cavity 2 can be related to the cavity mode a2 as

By substituting Eqs. (7) and (8) into Eqs. (5) and (6), we can obtain the expected value of the atomic operators, which takes the form

Here, we denote g1 ≡ |g1 |ei ϕ1 and g2 ≡ |g2 |ei ϕ2, Kjκj / 2 + i Δc,j (j = 1,2), Kaγa/2 + iΔa, and PKa (K1 K2 + J2) + K2 |g1|2 + K1 |g2|2.

In order to measure the squeezing, the field operator must be expressed in terms of quadrature components. For example, the output operator Sout should be divided into two quadrature components as[42]

Here θ is an adjustable phase, and we set θ = 0 below. XAmp and XPh are the amplitude and phase quadrature of the output field, respectively, which satisfy the commutation relation [XAmp,XPh] = i/2. According to the definition of the variance ⟨M,N⟩ = ⟨MN⟩ − ⟨M⟩ ⟨N⟩, i.e., ⟨Sout,Sout⟩ = ⟨Sout Sout ⟩ − ⟨Sout ⟩ ⟨Sout⟩, the corresponding quadrature variance of the output field Sout read as

The notation ⟨: , :⟩ represents the normally ordered variance, which originate from the definition ⟨: Xi,Xi :⟩ = ⟨Xi,Xi⟩ −1/4, i = Amp,Ph. If ⟨ :Xi,Xi:⟩ is negative, the output field exhibits quadrature squeezing.[42] In detail, ⟨:XAmp,XAmp:⟩ measures the squeezing of the amplitude quadrature, while ⟨:XPh,XPh:⟩ measures the squeezing of the phase quadrature.

The variances of the output operator Sout, i.e., ⟨Sout,Sout⟩ and , can be obtained from Eqs. (7)–(11),

It is clear from Eqs. (13)–(16) that the normally ordered variance of the output field is related to both the amplitude and phase of the coupling strengths g1 and g2. As a result, we can explore the phase-controlled squeezing effect in this system. Furthermore, in this work, we are mainly interested in the influence of the phase effect on the squeezing behaviors of the output field. Thus, in the following discussion, for simplicity but without loss of generality, we assume that the two cavities have the same frequency ωc,1 = ωc,2ωc and interact resonantly with the two-level atom ωa = ωc, i.e., Δc,1 = Δc,2 = Δa.

3. Results and discussion

In this section, we study the squeezing properties of the transmitted field and focus on discussing how the quadrature squeezing can be affected by the phase difference of the coupling strengths between the two-level atom and the two cavities. First, we discuss the case where the two cavities have the same decay rates κ1 = κ2 = κ and couple to the atom with the same coupling strength |g1| = |g2| = |g|. Then, the mismatch between the coupling strengths |g1| and |g2|, and the mismatch between the decay rates κ1 and κ2 will be considered later.

In Figs. 2(a)2(c), the variances of the transmitted field versus the detuning Δ/κ = (ωωc)/κ between the probe field and the cavity modes with different phase differences ϕ2ϕ1 are plotted. Here, we study how the variances of the amplitude and phase quadrature vary with the detuning, and show how the squeezing properties are modulated by the relative phase of the coupling strengths. It can be seen from Figs. 2(a)2(c) that there are several dips with variance below zero in the squeezing spectra. The number and location of the dips are strongly related to phase difference. If the atom couples to the two cavities in phase (ϕ2ϕ1 = 0), there is only a dip in the phase quadrature located at Δ = J, which means that only the phase quadrature can be squeezed, as shown in Fig. 2(a). If the atom couples to the two cavities out of phase (ϕ2ϕ1 = π), a similar squeezing property in the phase quadrature can be obtained, except that the location of the dip changes from Δ = J to Δ = −J, as shown in Fig. 2(b). However, If the atom couples to the two cavities with a relative phase π/2(ϕ2ϕ1 = π/2), the squeezing behaviors of the phase and amplitude quadrature are much different from that of the cases in phase or out of phase. There are two dips with equal variance in the phase quadrature located at Δ = − J and Δ = J, respectively. Furthermore, there is also a dip in the amplitude quadrature located at Δ = 0, which means that, in the present case, an obviously squeezing behavior also can be observed in the amplitude quadrature, as shown in Fig. 2(c).

Fig. 2. The variances of the amplitude and phase quadrature as a function of the detuning for different phase difference: (a) ϕ2ϕ1 = 0; (b) ϕ2ϕ1 = π; (c) ϕ2ϕ1 = 0.5π. All the red solid and blue dashed curves represent the amplitude and phase quadrature, respectively. The common parameters for this situation are |g1| = |g2| = 5κ, J = 2κ, ε = 0.5κ, γa = 0.1κ, and κ1 = κ2 = κ.

The physical mechanism of this phase-dependent quadrature squeezing can be explained in the normal mode picture.[43,44] The system contains two cavities interacting with an atom supports a pair of standing wave modes A and B, which are the superposition of the two cavity modes, i.e., with frequency ωc + J, and with frequency ωcJ. For convenience, we use the words high-frequency mode and low-frequency mode to describe the standing wave modes A and B, respectively. In this new normal mode representation, the Hamiltonian in Eq. (2) can be rewritten as

Here, Δϕ = ϕ2ϕ1. From Eq. (17), we can find that the effective coupling strengths between the atom and the standing wave modes A and B are given by and , respectively, and these effective coupling strengths can be modified by tuning the relative phase (Δ ϕ = ϕ2ϕ1). When |g1| = |g2| = |g|, the effective coupling strengths geff,A and geff,B become and , which can be as large as or as small as zero when the relative phase changes from zero to π. For the case ϕ2ϕ1 = 0, we can obtain and geff,B = 0, which means that the atom couples to the high-frequency mode A with a maximal effective coupling strength, while decouples completely from the low-frequency mode B. The coherent interaction between the atom and the normal mode A induces a strong nonlinearity in the system, which leads to the phase quadrature squeezing around the high-frequency ωc + J region (i.e., Δ = J). Vice versa, for the case ϕ2ϕ1 = π, the atom only couples to the normal mode B (geff,A = 0 and ), which result in a similar squeezing behavior around the low-frequency ωcJ region (i.e., Δ = − J). However, for the situation ϕ2ϕ1 = π/2, the atom couples to both the normal modes A and B. As a consequence, the phase quadrature squeezing can be observed simultaneously around the low-frequency and high-frequency region. Furthermore, the squeezing behavior in the amplitude quadrature can also be observed around the bare cavity resonant region Δ = 0 due to the joint contribution of the two normal modes A and B.

A remarkable feature of the quadrature squeezing of the transmitted field is that the squeezing behaviors can be transferred from the low-frequency mode to the high-frequency mode, and can be transferred from the phase quadrature to the amplitude quadrature by modulating the relative phase of the coupling strengths. To clearly show this, we first plot the variances of the phase and amplitude quadrature versus the phase difference with three special detuning Δ = ± J and Δ = 0 in Figs. 3(a)3(c), respectively. As the input field resonant with the high-frequency mode A (i.e., Δ = J), the squeezing behavior can be observed in the phase quadrature when the atom couples to the two cavities in phase ϕ2ϕ1 = 0. However, the degree of squeezing decreases with increasing ϕ2ϕ1. The squeezing behavior around the high-frequency mode A even disappear as the phase difference increases to π, as shown in Fig. 3(a). On the contrary, as the input field resonant with the low-frequency mode B (i.e., Δ = −J), no obvious squeezing behaviors around the low-frequency mode B can be obtained in both the amplitude and phase quadrature at first (ϕ2ϕ1 = 0). However, by increasing the phase difference to π, the phase quadrature can be squeezed and the degree of squeezing increases to the same value (about −0.2) as that obtained around the high-frequency mode A, as shown in Fig. 3(b). It is clear from Figs. 3(a) and 3(b) that the squeezing behaviors can be transferred completely from the high-frequency region to the low-frequency region with the phase difference ϕ2ϕ1 ranges from zero to π. Furthermore, figure 3(c) shows that, as the input field resonant with the bare cavity mode (i.e., Δ = 0), the squeezing behaviors can be transferred between the amplitude and phase quadrature by increasing the phase difference, and the squeezing behavior in the amplitude quadrature can only be obtained obviously with the phase difference around π/2.

Fig. 3. The variances of the amplitude and phase quadrature as a function of the phase difference for different detunings: (a) Δ = ωωc = J; (b) Δ = ω - ωc = − J; (c) Δ = ωωc = 0. All the red solid and blue dashed curves represent the amplitude and phase quadrature, respectively. The common parameters are the same as those in Fig. 2.

In order to further illuminate these features and study the effect of the relative phase on the quadrature squeezing in a more general case, the variances of the phase and amplitude quadrature versus both the detuning Δ/κ = (ωωc)/κ and the phase difference ϕ2ϕ1 are plotted in Figs. 4(a) and 4(b), respectively. We can find that, for a given detuning, the squeezing behaviors in the high and low-frequency regions are much different. The squeezing effect in the low-frequency region becomes more and more significant while that in the high-frequency region disappears gradually as the relative phase increasing, and vice versa. This shows the transition of the squeezing between the high and low-frequency regions, which is consistent with the results shown in Figs. 3(a) and 3(b). Furthermore, comparing Figs. 4(a) and 4(b), it is found that the squeezing behavior in the phase quadrature is also much different from that in the amplitude quadrature around the bare cavity resonant region (Δ/κ = (ωωc)/κ = 0). The variance of the phase or amplitude quadrature increases above zero or decreases blew zero alternatively with the change of the phase difference around π/2, which means that the squeezing between the phase quadrature and the amplitude quadrature can be transferred back and forth. Therefore, the squeezing transfer between the phase and amplitude quadrature or between the low and high-frequency regions can be achieved by simply modulating the relative phase of the coupling strengths between the atom and the two cavities.

Fig. 4. The variances of the phase quadrature (a) and the amplitude quadrature (b) as a function of both the phase difference and the detuning. The common parameters are the same as those in Fig. 2.

In the above discussion, we have assumed that the two cavities with the same decay rates κ1 = κ2, and coupled symmetrically to the two-level atom, i.e., |g1| = |g2|. We find that, under this conditions, the squeezing properties of the transmitted field depend strongly on the phase difference of the coupling strengths between the atom and the two cavities. However, it is clear from Eqs. (15) and (16) that the phase-dependent quadrature squeezing is also related to the intensity of the decay rates and the coupling constants. Below, we focus on analyzing the effect of decay mismatch κ1κ2 and coupling mismatch |g1| ≠ |g2| on the phase-dependent quadrature squeezing properties.

First, we study how the decay mismatch between the two cavities affects the quadrature squeezing by setting |g1| = |g2|. In Figs. 5(a)5(c), the variances of both the amplitude and phase quadrature are plotted as a function of the decay mismatch κ1/κ2 for three detunings, i.e., Δ = 0 and Δ = ± J. Here, we take the phase difference ϕ2ϕ1 = π/2 for example. Figure 5 shows clearly that the degree of squeezing depends strongly on the decay mismatch. As the decay rate ratio κ1/κ2 increases, both the variance of the phase quadrature located at Δ = ± J, and the variance of the amplitude quadrature located at Δ = 0 increase monotonically, as shown by the blue dashed lines in Fig. 5(a) and Fig. 5(c), and the red solid line in Fig. 5(b), respectively. Compare with Fig. 2(c), where κ1/κ2 = 1 (κ1 = κ2), the degree of squeezing increases when κ1/κ2 < 1 (κ1 < κ2). The reason is that, as the probe field is injected into cavity 1, an unexpected noise will be injected into the cavities at the same time, which weaken the quadrature squeezing of the transmitted field. However, as the decay rate of cavity 1 κ1 decreases, partial fluctuation in the input field can be effectively suppressed, which leads to the enhancement of the squeezing of the output field.[45] As a result, a better quadrature squeezing of the output field can be obtained with a proper choice of the decay mismatch κ1/κ2.

Fig. 5. The variances of the amplitude and phase quadrature as a function of the decay mismatch for different detunings: (a) Δ = ωωc = J; (b) Δ = ωωc = 0; (c) Δ = ωωc = − J. All the red solid and blue dashed curves represent the amplitude and phase quadrature, respectively. The common parameters for this situation are |g1| = |g2| = 5κ2, J = 2κ2, ε = 0.5κ2, and γa = 0.1κ2.

Then, we proceed to study how the coupling mismatch between the two cavities and the atom affects the quadrature squeezing by setting κ1 = κ2. In Figs. 6(a) and 6(b), the quadrature squeezing is plotted as a function of the coupling mismatch |g2|/|g1| with the phase difference ϕ2ϕ1 = 0 for two detuning, i.e., Δ = ± J. For the case ωωc = − J, the variances of the phase and amplitude quadrature are approach to zero when |g2|/|g1| = 1 (|g1| = |g2|), which means that no obvious squeezing effect can be observed around Δ = − J. This is consistent with the result shown in Fig. 2(a). However, when |g2|/|g1| ≠ 1 (|g1| ≠ |g2|), the variances of both the phase and amplitude quadrature can become much smaller than zero, which means that a better quadrature squeezing can be obtained due to the asymmetric interaction between the atom and the two cavities. Furthermore, the squeezing behaviors of the phase and amplitude quadrature can be transferred into each other by modulating the coupling mismatch, i.e., the squeezing effect of the phase quadrature can be transferred into the amplitude quadrature by increasing |g2|/|g1| at first, and then transferred back to the phase quadrature by further increasing the coupling mismatch, as shown in Fig. 6(a). While for the case ωωc = J, the squeezing behaviors can only be observed in the phase quadrature, and the degree of squeezing is not very sensitive to the coupling mismatch |g2|/|g1|, as shown in Fig. 6(b).

Fig. 6. The variances of the amplitude and phase quadrature as a function of the coupling mismatch for two detunings: (a) Δ = ωωc = − J; (b) Δ = ωωc = J. All the red solid and blue dashed curves represent the amplitude and phase quadrature, respectively. The common parameters for this situation are J = 2κ, ε = 0.5κ, γa = 0.1κ, and κ1 = κ2 = κ.

Compared to Fig. 2(a), we can find that the squeezing behaviors of the transmitted field become much different due to the coupling mismatch. This difference can be further explaining in the normal mode picture. As mentioned above, for the case |g1| = |g2| = |g|, the effective coupling strengths between the atom and the standing wave modes A and B, i.e., and , can reach to their maximum value and minimum value geff,B = 0, respectively, as Δϕ = ϕ2ϕ1 = 0. This result in an obvious quadrature squeezing around the high-frequency mode A (Δ = J), but no quadrature squeezing can be observed around the low-frequency mode B (Δ = −J), as shown in Fig. 2(a). However, for the situation |g1| ≠ |g2|, the expressions of the effective coupling strengths should be rewritten as and , respectively. It is clear that, in the present case, the effective coupling strength geff,B becomes nonzero even though the atom couples to the two cavities in phase (Δϕ = 0), which leads to the quadrature squeezing around the low-frequency mode B (Δ = −J), as shown in Fig. 6(a).

Furthermore, it is worth noting that the effect of the coupling mismatch on squeezing with other phase difference can be studied in the same way as those for Δϕ = 0. Actually, the quadrature squeezing vary with |g2|/|g1| for the phase difference Δϕ = π also have been investigated. It is confirmed that a similar squeezing behavior can be obtained and the results are the mirror images of those in Fig. 6. Namely, as the two cavities couple to the atom out of phase Δϕ = π, the squeezing behaviors and the degree of squeezing located at the high-frequency mode A (Δ = J) are very sensitive to |g2|/|g1|. Obviously, a better quadrature squeezing effect can also be achieved by adjusting the coupling mismatch between |g1| and |g2|.

Before summary, let us provide a brief discussion on the experimental feasibility of our scheme. The interaction between a two-level atom and two directly coupled cavities can be realized by coupling a quantum dot to a semiconductor bimodal microcavity,[46,47] which supports two degenerate counterpropagating cavity modes (a clockwise propagating mode and a counterclockwise propagating mode). The coupling between these two counterpropagating modes can be realized due to the surface scattering of a nanoparticle, and the mode coupling strength can be controlled by tuning the radius of the nanoparticle or the number of the nanoparticles.[48] Furthermore, in our model, a key element is that the atom couples to the two cavity modes with different relative phase. It has been demonstrated that, in the bimodal microcavity, the field distributions of the clockwise and counterclockwise propagating modes have an azimuthal dependence.[43,49] Thus, the coupling strengths and their relative phase between the two-level atom and the two cavities can be modulated experimentally by varying the position of the atom relative to the azimuthal distribution of the cavity modes.

4. Conclusion

In this paper, we have investigated the quadrature squeezing of the output field from a cavity QED system, which contains two coupled cavities and a two-level atom. We mainly focus on studying the influence of the phase effect of the coupling strengths on the squeezing. It is found that the squeezing behaviors can be transferred between the phase quadrature and the amplitude quadrature, and can be transferred between the low-frequency region and high-frequency region by tuning the phase difference of the coupling strengths between the atom and the two cavities. Besides, the decay mismatch between the two cavities and the coupling mismatch between the atom and the two cavities can be used to generate better quadrature squeezing. The degree of squeezing of the output field can be remarkably enhanced with the help of these mismatch, which is useful for the application in highly precision measurements and atomic spectroscopy.

Reference
[1] Nielsen M A Chuang I L 2011 Quantum Computation and Quantum Information Cambridge Cambridge University Press
[2] Houck A A Türeci H E Koch J 2012 Nat. Phys. 8 292
[3] Reiserer A Rempe G 2015 Rev. Mod. Phys. 87 1379
[4] Gisin N Thew R 2007 Nat. Photon. 1 165
[5] Kok P Munro W J Nemoto K Ralph T C Dowling J P Milburn G J 2007 Rev. Mod. Phys. 79 135
[6] Carmichael H J 1985 Phys. Rev. Lett. 55 2790
[7] Pan J W Chen Z B Lu C Y Weinfurter H Zeilinger A Zukowski M 2012 Rev. Mod. Phys. 84 777
[8] Andersen U L Gehring T Marquardt C Leuchs G 2016 Phys. Scr. 91 053001
[9] Loudon R Knight P L 1987 J. Mod. Opt. 34 709
[10] Yang S Zhang S Zou X Bi S Lin X 2013 Phys. Rev. 87 024302
[11] Abadie J Abbott B P Abbott R et al. 2011 Nat. Phys. 7 962
[12] Szigeti S S Tonekaboni B Lau W Y S Hood S N Haine S A 2014 Phys. Rev. 90 063630
[13] Lvovsky A I Raymer M G 2009 Rev. Mod. Phys. 81 299
[14] Slusher R E Hollberg L W Yurke B Mertz J C Valley J F 1985 Phys. Rev. Lett. 55 2409
[15] Wu L A Kimble H J Hall J L Wu H 1986 Phys. Rev. Lett. 57 2520
[16] Vahlbruch H Mehmet M Danzmann K Schnabel R 2016 Phys. Rev. Lett. 117 110801
[17] Li S Li H Z Xu J P Zhu C J Yang Y P 2019 Acta Phys. Sin. 68 174202 in Chinese
[18] Rashid M Tufarelli T Bateman J Vovrosh J Hempston D Kim M S Ulbricht H 2016 Phys. Rev. Lett. 117 273601
[19] Agarwal G S Huang S M 2016 Phys. Rev. 93 043844
[20] Német N Parkins S 2016 Phys. Rev. 94 023809
[21] Adamyan H H Bergou J A Gevorgyan N T Kryuchkyan G Y 2015 Phys. Rev. 92 053818
[22] Meystre P Zubairy M S 1982 Phys. Lett. 89 390
[23] Wang C Q Zou J Zhang Z M 2016 Chin. Phys. Lett. 33 024202
[24] Li H Xu J Wang D W Xia X Yang Y Zhu S 2017 Phys. Rev. 96 013832
[25] Tian J F Zuo G H Zhang Y C Li G Zhang P F Zhang T C 2017 Chin. Phys. 26 124206
[26] Li Z Ma S L Li F L 2015 Phys. Rev. 92 023856
[27] Liu F F Fang M F Xu X 2019 Chin. Phys. 28 060304
[28] Chuang Y L Lee R K Yu I A 2017 Phys. Rev. 96 053818
[29] Ourjoumtsev A Kubanek A Koch M Sames C Pinkse P W H Rempe G Murr K 2011 Nature 474 623
[30] Schulte C H H Hansom J Jones A E Matthiesen C Gall C L Atature M 2015 Nature 525 222
[31] Kim N C Li J B Yang Z J Hao Z H Wang Q Q 2010 Appl. Phys. Lett. 97 061110
[32] Xiao Y F Li M Liu Y C Li Y Sun X Gong Q 2010 Phys. Rev. 82 065804
[33] Cheng M T Xu J Agarwal G S 2017 Phys. Rev. 95 053807
[34] Chen G Y Lambert N Chou C H Chen Y N Nori F 2011 Phys. Rev. 84 045310
[35] Yang D C Cheng M T Ma X S Xu J Zhu C Huang X S 2018 Phys. Rev. 98 063809
[36] Yuan L Xu S Fan S 2015 Opt. Lett. 40 5140
[37] Pleinert M O Zanthier J V Agarwal G S 2018 Phys. Rev. 97 023831
[38] Wang C Liu Y L Wu R Liu Y X 2017 Phys. Rev. 96 013818
[39] Qu Y Li J Wu Y 2019 Phys. Rev. 99 043823
[40] Gardiner C W Collett M J 1985 Phys. Rev. 31 3761
[41] Auffèves-Garnier A Simon C Gérard J M Poizat J P 2007 Phys. Rev. 75 053823
[42] Collett M J Gardiner C W 1984 Phys. Rev. 30 1386
[43] Srinivasan K Painter O 2007 Phys. Rev. 75 023814
[44] Majumdar A Bajcsy M Rundquist A Vuckovic J 2012 Phys. Rev. Lett. 108 183601
[45] Scully M O Zubairy M S 1997 Quantum Optics Cambridge Cambridge University Press
[46] Srinivasan K Painter O 2007 Nature 450 862
[47] Min B Ostby E Sorger V Ulin-Avila E Yang L Zhang X Vahala K 2009 Nature 457 455
[48] He L Ozdemir S K Zhu J Monifi F Yilmaz H Yang L 2013 New J. Phys. 15 073030
[49] Aoki T Dayan B Wilcut E Bowen W P Parkins A S Kippenberg T J Vahala K J Kimble H J 2006 Nature 443 671